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The  results  of  a  series  of  numerical  simulations  are  presented  of  experiments  con¬ 
ducted  at  the  ONERA  Chalais-Meudon  Research  Center  and  at  the  CALSPAN-University 
at  Buffalo  Research  Center  on  shock-shock  interactions.  The  flowfield  characteristics  are 
described  with  the  aid  of  the  numerical  predictions  and  the  computed  values  of  surface 
pressure  and  heat  transfer  are  compared  with  the  experimental  measurements  for  pur¬ 
poses  of  verification  and  validation.  A  preliminary  grid  refinement  study  is  conducted  in 
some  test  conditions  to  investigate  certain  discrepancies  that  have  been  found  between 
the  numerical  and  the  experimental  data. 


Nomenclature 

Re  Reynolds  number  per  meter,  m_1 

M  Mach  number 

p  Density,  kg/m 3 

p  Pressure,  Pa 

T  Temperature,  K 

Rcyi  Cylinder  radius,  reference  lenght,  m 

x,  y  Cartesian  body  axes,  m 

Subscripts 

oo  Freestream  conditions 

w  Wall  conditions 

Superscripts 

o  Stagnation  conditions 

Introduction 

HE  correct  evaluation  of  the  heat  load  produced 
by  shock/shock  interaction  represents  an  impor¬ 
tant  and  critical  problem  in  the  design  of  hypersonic 
vehicles.  It  is  well  known,  in  fact,  that  shock/ shock 
interference  produces  very  high  levels  of  heat  transfer 
and  surface  pressure  in  the  regions  of  the  aircraft  that 
are  close  to  the  impinging  point.  In  air-breathing  hy¬ 
personic  vehicles,  the  phenomenon  is  typically  located 
at  the  intake  cowl  lip,  where  the  shocks  produced  by 
the  compression  ramps  intersects  the  bow  shock  ahead 
of  the  cowl,  but  it  may  affect  in  a  similar  manner  the 
wing  or  the  fin  leading  edge  also. 

In  a  fundamental  work  for  the  study  of  shock-on- 
shock  interactions,  Edney1  classified  six  types  of  flow 
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patterns,  which  depend  on  the  relative  position  of  the 
incoming  shock  with  respect  to  the  bow  shock.  Be¬ 
tween  them,  the  so-called  type  III  and  type  IV  interac¬ 
tion  are  the  most  critical,  since  they  are  characterized 
by  the  largest  values  of  wall  pressure  and  heat  trans¬ 
fer.  In  the  Edney  type  III  interaction,  pressure  and 
heat  flux  peaks  at  the  wall  are  related  to  the  impact 
of  the  shear  layer  produced  by  the  shock  interference, 
while  in  the  Edney  type  IV  interaction,  even  larger 
loads  are  due  to  the  impingement  of  a  supersonic  jet 
on  the  body  surface. 

In  the  last  years,  some  experimental  campaigns 
have  been  conducted  with  the  aim  of  better  under¬ 
standing  the  phenomenon  and  of  providing  data  for 
code  validation.  In  fact,  even  though  the  Compu¬ 
tational  Fluid  Dynamics  is  widely  used  today  in  the 
prediction  of  hypersonic  flows  for  design  purposes, 
it  is  still  uncertain  whether  it  can  be  considered  as 
reliable  for  simulating  any  kind  of  flow,  expecially 
when  they  present  complicated  features  and  different 
scales.  The  two  most  recent  sets  of  experiments  on 
shock/shock  interaction  were  conducted  at  ONERA 
Chalais-Meudon  (France)  in  the  wind  tunnel  R5Ch2 
and  at  the  Calspan-University  at  Buffalo  Research 
Center  (CUBRC)  (U.S.A.),  in  the  Large  Energy  Na¬ 
tional  Shock  (LENS)  tunnel.3  Subsequently,  various 
numerical  simulations  have  been  carried  out  for  re¬ 
producing  the  experimental  data,  in  particular  those 
obtained  by  ONERA. 4-6  In  addition,  both  the  exper¬ 
iments  have  been  selected  for  validating  CFD  capa¬ 
bility  in  the  framework  of  the  activities  of  the  NATO 
Research  Technology  Organization  (RTO)  Advanced 
Vehicle  Technology  (AVT)  Working  Group  10,  Sub¬ 
group  3  ”CFD  Validation  for  Hypersonic  Flight”.7 

In  this  paper,  the  results  of  numerical  predictions  of 
both  the  French  and  the  American  experiments  will  be 
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presented  and  a  validation  study  will  be  done. 

MODEL  SETUP  AND  TEST 
CONDITIONS. 

The  ONERA  experiments.2 

The  experiments  carried  out  at  ONERA  in  R5Ch 
blow-down  wind  tunnel  were  characterized  by  the 
nominal  stagnation  and  freestream  conditions  listed 
in  table  1.  The  experimental  set-up  that  produced 
the  Edney  type  IV  interaction  is  shown  in  figure  1. 
The  shock  generator,  a  prism  whose  cross  section  is 
an  isosceles  triangle  with  the  base  100  mm  long  and 
leading  edge  angle  of  10°,  has  a  spanwise  dimension  of 
100  mm.  In  the  present  experiments,  it  was  rotated  of 
10°  with  respect  to  the  freestream,  so  that  a  wedge/flat 
plate  configuration  with  20°  wedge  angle  was  obtained. 
The  cylinder,  which  was  placed  perpendicular  to  the 
freestream,  has  a  radius  of  8  mm.  Since  its  span  is  100 
mm  wide,  a  two-dimensional  flow  could  be  assumed 
at  the  middle  of  the  model,  where  the  measurement 
equipment  was  mounted.  The  relative  locations  of  the 
cylinder  and  the  wedge  could  be  varied  in  order  to 
obtain  different  types  of  shock-shock  interactions.  In 
particular,  for  the  Edney  type  IV  interference,  the  dis¬ 
tance  between  the  leading  edge  of  the  shock  generator 
and  the  cylinder  center  was  110  mm  along  the  x-axis 
and  53  mm  along  the  y-axis.  Two  models  of  the  cylin¬ 
der  were  built  and  equipped  with  pressure  taps  and 
thermocouples  to  measure  surface  pressure  and  heat 
transfer  respectively.  In  addition,  a  Dual-line  Coher¬ 
ent  Anti  Stokes  Scattering  (DL-CARS)  technique  was 

Table  1  Experimental  conditions  at  ONERA. 


J^O 

1050 

K 

p° 

=  250000 

Pa 

T 

- L  oo 

=  52.5 

K 

Poo 

=  5.9 

Pa 

Moo 

=  9.95 

Bp 

J  I'C'OO 

=  1.66  -105 

m_1 

Rcyl 

=  0.008 

m 

y* 


used  to  perform  instantaneous  measurements  of  tem¬ 
perature  and  density  inside  the  flowfield. 

The  CALSPAN  experiments.3 

Experimental  studies  were  conducted  at  CALSPAN 
in  the  Veridian  48-inch  shock  tunnel  to  obtain  detailed 
heat  transfer  and  pressure  measurements  in  regions  of 
shock/shock  interaction  in  laminar  and  turbulent  low- 
density  flows  at  Mach  numbers  from  10  to  16.  Seven 
runs  were  conducted  in  the  laminar  regime,  two  at 
Mach  15.7  using  a  cylinder  with  radius  of  0.375  inches 
and  five  at  Mach  14  using  a  cylinder  with  radius  of  1.5 
inches.  Other  seven  runs  were  in  the  turbulent  regime 
at  Mach  10.8,  using  a  cylinder  with  radius  of  1.5  inches. 
Miniature  high  frequency  instrumentation  was  used  to 
resolve  the  large  heat  transfer  gradients  and  the  flow 
unsteadiness  observed  in  these  studies.  In  order  to 
minimize  errors  in  the  heat  transfer  measurements, 
which  may  arise  from  lateral  heat  conduction  effects 
inside  the  model  in  case  of  large  peak  heating  values, 
it  was  preferred  to  use  models  with  non-conducting 
surfaces.  Then,  in  those  experiments  where  the  result¬ 
ing  wall  temperature  rises  were  important,  the  surface 
temperature  distribution  was  tabulated  together  with 
the  heat  transfer  rate,  so  that  it  could  be  used  in  future 
numerical  predictions  as  a  boundary  condition.3  The 
data  obtained  in  the  experimental  studies  have  been 
incorporated  in  a  database  called  CUBDAT,8  which 
can  be  accessed  by  researchers  involved  in  numerical 
prediction  studies.  In  this  paper,  numerical  results 
related  to  the  conditions  of  runs  38  and  43  of  the  CUB¬ 
DAT  database  will  be  presented. 

Table  2  Experimental  conditions  at  CALSPAN 
for  RUN  38  and  RUN  43 


RUN  38 

RUN  43 

J^O 

= 

3211.1 

3355.6 

K 

p° 

= 

3212957 

2992325 

Pa 

T 

-*■  oo 

= 

100.0 

110. 

K 

Poo 

= 

3.66 

3.51 

Pa 

Poo 

= 

1.273-10-4 

1.108-10-4 

kg/m3 

Moo 

= 

14.15 

13.94 

RCoo 

= 

52165 

42618 

m-1 

T 

-*■  w 

= 

296.67 

300. 

K 

Rcyl 

= 

0.0381 

0.0381 

m 

yi 


Fig.  1  Model  used  in  the  ONERA  experiments  for  Fig.  2  Model  used  in  the  CALSPAN  experiments 
the  EDNEY  type  IV  interaction  for  the  EDNEY  type  IV  interaction.  RUN  38. 


9-2 


RTO-TR-AVT-007-V3 


Bow  shock  Edney-type  I 


a)  Edney-type  I  interaction. 


Fig.  3  Sketches  of  the  Edney-type  I-II-III  and  IV  interactions. 


NUMERICAL  METHOD 

The  computational  results  that  will  be  presented  in 
this  paper  have  been  obtained  using  a  well  established 
numerical  technique,  whose  details,  which  will  be  not 
described  here,  can  be  found,  for  instance,  in  refer¬ 
ences,9.10  In  synthesis,  the  Navier-Stokes  equations 
for  two-dimensional  flows  are  solved.  The  numeri¬ 
cal  method  adopts  a  finite  volume  discretization  of 
the  physical  domain,  an  explicit  integration  scheme,  a 
flux  difference  splitting  technique11  for  the  evaluation 
of  the  convective  fluxes  and  a  standard  cell  centered 
method  for  the  diffusive  fluxes.  A  multi-block  do¬ 
main  decomposition  method  is  used,  that  allows  for 
an  efficient  distribution  of  grid  points  across  the  com¬ 
putational  domain.  Finally,  a  nominal  second  order 
accuracy  in  space  and  time  is  reached  using  an  Essen¬ 
tially  Non-Oscillatory  scheme.12 

FLOWFIELD  FEATURES 

Depending  on  the  strength  of  both  the  impinging 
and  the  bow  shocks  and  on  their  relative  positions, 
various  shock  interference  patterns  can  occur.  In  his 
fundamental  work,  Edney1  classified  six  types  of  in¬ 
teractions.  The  first  four  interaction  patterns  are 
sketched  in  figure  3  for  comparison  with  the  numer¬ 
ical  results  that  will  be  shown  in  the  next  lines  and 
the  reader  is  encouraged  to  consult  the  Edney’s  report 
for  details. 

In  the  following,  the  results  of  several  numerical 
tests  performed  by  the  author  changing  the  incoming 
shock  position  will  be  shown  in  order  to  demonstrate 


the  effect  of  the  various  interference  patterns  on  the 
flowfield  and  on  the  wall  pressure  and  heat  transfer 
levels.  The  freestream  conditions  and  the  incoming 
shock  strength  are  those  of  the  CALSPAN’s  run  38. 
The  mesh  used  in  such  computations  is  composed  of 
150  points  in  the  normal  direction,  with  grid  stretch¬ 
ing  close  to  the  wall.  In  the  tangential  direction,  the 
grid  spacing  is  1  degree  in  the  range  [—40°  :  +40°]. 
Then,  the  tangential  grid  spacing  increases  linearly  up 
to  a  maximum  value  of  1.5  degrees  at  the  upper  and 
lower  boundaries.  The  Mach  number  contours  and  the 
streamlines  patterns  are  shown  in  figure  4  for  twelve 
different  shock  positions.  In  the  case  of  the  first  shock 
position  (#1)  an  Edney-type  IV  interaction  occurs, 
but  the  resulting  supersonic  jet  does  not  impinge  on 
the  cylinder  wall.  The  stagnation  point  is  very  close 
to  the  symmetry  axis,  but,  as  it  will  be  shown  later, 
the  pressure  and  heat  flux  peak  values  are  higher  with 
respect  to  the  undisturbed  cylinder,  since  the  stag¬ 
nation  streamline  comes  from  the  high  pressure  side 
of  the  impinging  shock  wave.  Lowering  the  incom¬ 
ing  shock  position  (#2  to  #8),  the  supersonic  jet  is 
directed  towards  the  body  surface  and  the  flow  ex¬ 
periences  a  final  strong  compression  across  a  normal 
shock  at  a  short  distance  from  the  wall.  Lowering  the 
shock  a  little  more  (#9),  the  slip  surface  grazes  the 
wall  and  the  normal  shock  in  front  of  the  body  dis¬ 
appears.  This  is  an  Edney-type  III  interaction.  If  we 
continue  to  move  the  incoming  shock  downwards  (#10 
and  #11),  we  would  expect  to  see  an  Edney-type  II  in¬ 
teraction.  In  these  particular  two  cases,  however,  the 
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a)  Incoming  shock  position  #1.  b)  Incoming  shock  position  #2.  c)  Incoming  shock  position  #3. 


d)  Incoming  shock  position  #4.  e)  Incoming  shock  position  #5.  f)  Incoming  shock  position  #6. 


g)  Incoming  shock  position  #7.  h)  Incoming  shock  position  #8.  i)  Incoming  shock  position  #9. 


j)  Incoming  shock  position  #10.  k)  Incoming  shock  position  #11.  1)  Incoming  shock  position  #12. 

Fig.  4  Effect  of  the  impinging  shock  position.  Conditions  of  CALSPaN’s  RUN  43.  Mach  number 
contours. 
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a)  Pressure  at  the  wall. 


b)  Heat  flux  at  the  wall. 

Fig.  5  Pressure  and  heat  flux  at  the  wall  for  dif¬ 
ferent  incoming  shock  positions.  Conditions  of  the 
CALSPAN’s  run  43. 

reflected  shock  R  interacts  with  the  subsonic  pocket 
on  the  cylinder  nose  and  the  Edney-type  II  interac¬ 
tion  is  not  possible.  Finally,  an  even  lower  impinging 
shock  (#12)  produces  and  Edney-type  I  interaction. 
Pressure  and  heat  transfer  levels  at  the  wall  for  all  the 
twelve  situations  are  shown  in  figure  5a  and  5b.  The 
Edney-type  IV  interaction  clearly  produces  the  high¬ 
est  peak  levels,  in  particular  when  the  jet  is  directed 
towards  the  cylinder  and  the  normal  shock  is  present 
in  front  of  the  body.  The  highest  pressure  levels  are 
related  to  condition  #3  and  #8,  which  correspond  to 
a  situation  when  the  stagnation  streamline  crosses  the 
final  normal  shock  respectively  at  its  upper  and  lower 
ends,  where  it  is  less  intense. 

RESULTS 

Numerical  simulation  of  the  ONERA  experiments 

The  experiment  conducted  at  the  ONERA  concern¬ 
ing  the  Edney  type  IV  interaction  has  been  replicated 
with  a  numerical  simulation  using  the  freestream  con¬ 
ditions  listed  in  table  3.  In  predicting  this  experiment, 


Table  3  Numerical  freestream  conditions  for  sim¬ 
ulating  the  ONERA  experiment. 


Poo 

5.9 

Pa 

Moo 

=  9.95 

Uoo 

=  1450 

m/s 

Too 

=  52.5 

K 

poo 

=  3.910  TO-4 

kg/m3 

poo 

=  3.405  T0-® 

Pa-s 

RCqo 

=  1.66  TO5 

Rcyl 

=  0.008 

m 

a  particular  care  must  be  employed  in  correctly  sim¬ 
ulating  the  flow  about  the  shock  generator  also.  In 
fact,  due  to  the  model  set-up  used  (figure  1),  an  in¬ 
tense  expansion  is  generated  at  the  corner  between 
the  upstream  part  of  the  shock  generator,  which  is 
a  wedge,  and  the  second  part,  which  is  a  flat  plate. 
The  expansion  fan  bends  the  oblique  shock  wave  and 
makes  the  flow  in  high-pressure  side  of  it  not  uniform. 
In  addition  to  this,  the  smallness  of  the  cylinder  ra¬ 
dius  with  respect  to  the  shock  generator  characteristic 
length  implies  that  a  slight  variation  of  the  inclination 
(or  position)  of  the  impinging  shock  may  produce  a  no¬ 
ticeable  displacement  in  the  position  of  the  stagnation 
point.  Therefore,  for  validation  purposes,  it  is  neces¬ 
sary  to  numerically  reproduce  with  great  accuracy  the 
flow  about  the  shock  generator. 

Computations  have  been  performed  using  one  grid 
for  the  upper  part  of  the  shock  generator,  one  grid 
for  the  lower  part,  one  grid  for  the  stream  past  the 
shock  generator  and  finally  a  grid  composed  of  300x351 
points  in  the  normal  and  tangential  directions,  respec¬ 
tively,  for  the  cylinder.  All  the  grids  are  stretched  in 
the  normal  direction  close  to  the  solid  walls.  The  mesh 
around  the  cylinder  is  refined  with  stretching  also  in 
the  region  of  the  supersonic  jet.  The  external  bound¬ 
aries  of  the  four  grids  used  are  visible  in  figure  6. 

In  order  to  estimate  the  scale  of  the  grid  cells  in 
relation  with  flow  conditions,  the  dimensionless  dis¬ 
tance  from  the  wall  y+  of  the  first  cell  centers  has 
been  computed.  The  plot,  not  included  here  for  lack 
of  space,  shows  that  first  cells  centers  lay  at  distances 
y+  <  0.2  for  most  part  of  the  cylinder  walls,  and  reach 
the  maximum  value  of  0.74  in  correspondence  with 
the  stagnation  point.  Thanks  to  this  analysis  and  to 
preliminary  investigations  with  coarser  grids,  we  are 
confident  that  the  grid  convergence  has  been  reached. 

In  this  particular  test  case,  the  intersection  between 
the  incident  oblique  shock  generated  by  the  wedge  and 
the  detached  shock  which  forms  ahead  of  the  cylinder 
takes  place  where  the  bow  shock  is  strong  (subsonic 
flow  in  the  high  pressure  side).  The  interaction  belongs 
to  the  Edney-type  IV  family1  and  is  characterized  by 
the  presence  of  a  supersonic  jet  surrounded  by  subsonic 
flow.  In  principle,  a  repeated  pattern  of  expansions 
and  compressions  develops  inside  the  jet  channel,  but 
in  this  case,  due  to  the  vicinity  of  the  jet  to  the  cylinder 
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a)  LogiO  of  pressure  contours. 
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b)  Mach  number  contours. 

Fig.  6  Simulation  of  the  complete  flowfield  for  the 
ONERA  experiment. 

wall,  we  can  distinguish  only  a  few  of  them.  In  partic¬ 
ular,  from  the  analysis  of  our  results,  we  interpret  the 
flowfield  inside  the  jet  as  sketched  in  figure  7b).  The 
shock  wave  at  the  ”  entrance”  the  supersonic  channel  is 
irregularly  reflected  from  the  upper  slip  line  and  forms 
a  A-shock.  In  correspondence  with  the  reflection  point, 
an  expansion  fan  forms  to  permit  to  part  of  the  stream¬ 
lines  to  flow  slightly  upwards.  At  the  opposite  side  of 
the  jet,  the  second  leg  of  the  lambda- shock  interacts 
with  the  lower  slip  line  and  is  reflected  as  an  expansion 
fan,  which  deviates  part  of  the  jet  slightly  downwards. 
Approaching  the  body,  the  supersonic  stream  which 
flows  inside  the  jet  is  suddenly  decelerated  through 
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Fig.  7  Investigation  of  the  wave  pattern  inside  the 
supersonic  jet. 

a  detached  shock.  The  jet  is  finally  divided  in  two 
streams  which  flow  upward  and  downward  the  stag¬ 
nation  point.  The  previous  interpretation  arise  from 
the  study  of  figure  7,  where  logio  of  pressure  plots  are 
shown  in  part  a),  with  the  paths  of  two  streamlines 
superimposed.  In  figures  7c)  and  7d),  the  pressure  con¬ 
ditions  along  these  streamlines  are  displayed.  It  can  be 
noticed  that,  after  the  pressure  rise  due  to  the  oblique 
Mach  shock,  a  second  increase  in  pressure,  due  to  the 
shock  at  the  ’’entrance”  of  the  channel,  is  present. 
This  second  pressure  rise  is  stronger  and  unique  for 
the  upper  streamline,  which  crosses  the  A-shock  prac¬ 
tically  at  the  triple  point,  while  it  is  less  intense,  and 
followed  by  a  second  compression  for  the  lower  stream¬ 
line,  which  crosses  the  two  legs  of  the  A-shock.  Then, 
both  streamlines  experience  a  pressure  drop  due  to 
the  expansions  fans  and  finally  a  strong  compression 
through  the  detached  shock  in  front  of  the  body.  This 
last  pressure  rise  is  smoother  for  the  upper  streamline 
than  for  the  lower  one,  because  the  former  crosses  the 
detached  shock  at  its  margins,  where  it  is  just  starting 
to  appear,  created  by  converging  compression  waves. 

A  larger  scale  view  of  the  interaction  is  shown  in 
figure  8.  Notice  that  the  position  of  the  cylinder  bow 
shock  is  completely  different  with  respect  to  the  undis¬ 
turbed  case,  being  the  shock  layer  thickness  of  the 
upper  portion  more  than  doubled. 

Computed  pressure  and  heat  flux  distributions  at 
the  wall  are  shown  in  figure  9,  where  they  are  also 
compared  with  results  of  experiments  conducted  in  the 
ONERA  R5Ch  low  density  wind  tunnel.4  The  com- 


a)  Logio  of  pressure  con-  b)  Mach  number  con¬ 
tours.  tours. 

Fig.  8  The  interaction  region. 
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theta  [deg] 

b)  Surface  heat  transfer  distribution. 


Fig.  9  Pressure  and  heat  flux  distribution  at  the 
wall.  Experimental  values  versus  numerical  predic¬ 
tions. 

parison  appears  to  be  rather  good.  Concerning  the 
level  of  the  pressure  peak,  we  explain  the  higher  value 
displayed  by  the  numerical  simulation  with  the  bad 
spatial  resolution  of  the  experiments,  due  to  the  fact 
that  in  the  latter  the  diameter  of  the  pressure  holes  is 
rather  large  (1.5  mm).13  The  situation  is  improved  for 
the  heat  fluxes,  since  the  diameter  of  the  thermocou¬ 
ples  is  smaller  (0.4  mm). 

In  order  to  simulating  the  effect  of  the  dimension 
of  the  pressure  taps  and  of  the  thermocouples  on  the 
results,  the  numerical  data  have  been  averaged  out 


theta  [deg] 

b)  Surface  heat  transfer  distribution. 

Fig.  10  Pressure  and  heat  flux  distribution  at  the 
wall.  Experimental  values  versus  averaged  and  not 
averaged  numerical  predictions. 

around  every  experimental  measurement  point  con¬ 
sidering  a  pressure  tap  with  a  1.5  mm  diameter  and 
a  thermocouple  with  a  0.4  mm  diameter.  For  the 
pressure  taps,  this  results  in  averaging  out  across  a 
circumferential  angle  of  ±5.37°  with  respect  to  the 
nominal  measurement  point.  A  weighted  average  has 
been  done  to  keep  into  account  the  stretching  of  the 
grid  in  the  tangential  direction.  The  results  obtained 
for  the  surface  pressure  distribution  are  very  interest¬ 
ing  (figure  10).  In  fact,  one  can  see  that  the  peak 
level  is  the  same  when  comparing  experiments  and 
(averaged)  numerical  predictions.  The  improvement 
is  definitely  less  important  when  observing  the  surface 
heat  flux  distribution,  but  in  this  case  the  experimen¬ 
tal  data  are  few,  and  it  is  impossible  to  decide  whether 
the  largest  experimental  is  the  real  peak  or  not. 

Numerical  simulation  of  the  CALSPAN 
experiments 

In  this  section,  numerical  results  related  to  RUN  38 
and  RUN  43  will  be  shown  and  compared  with  the  ex¬ 
perimental  data.  The  conditions  for  the  two  test  cases 
are  listed  in  table  2.  For  both  cases  the  correct  incom¬ 
ing  shock  position  to  match  with  the  experiments  was 
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Shock-shock  interactions  -  RUN  38  -  CALSPAN 


a)  Mach  number  contours. 


b)  Surface  heat  transfer  distribution. 

Fig.  11  Mach  number  and  logio  of  pressure  for 
RUN  38. 

obtained  moving  the  shock  until  the  wall  pressure  and 
heat  flux  peaks  were  in  the  same  tangential  position 
of  the  experimental  data. 

Concerning  RUN  38,  the  computations  were  per¬ 
formed  using  two  different  grids.  The  coarser  grid  is 
composed  of  150  points  in  the  normal  direction,  with 
a  constant  spacing  of  1  degree  in  the  tangential  di¬ 
rection.  The  finer  mesh  has  300  points  normal  to  the 
wall  and  the  constant  tangential  spacing  is  0.5  degrees. 
The  resulting  configuration  is  an  Edney-type  IV  inter¬ 
action,  as  it  is  evident  from  figures  11a  and  lib.  The 
effect  of  the  grid  refinement  on  the  wall  pressure  and 
heat  flux  is  shown  in  figures  12a  and  12b,  where  also 
the  experimental  data  are  plotted.  As  it  typically  hap¬ 
pens  in  the  numerical  simulation  of  this  kind  of  flows, 
refining  the  mesh  increases  the  peak  levels.  The  reason 
for  such  a  behaviour  is  not  clear,  but  again  it  could  be 


Angle  [deg.] 

a)  Surface  pressure  distribution. 


Shock-shock  interactions  -  RUN  38  -  CALSPAN 


Fig.  12  Pressure  and  heat  flux  distribution  at  the 
wall  for  RUN  38.  Experimental  values  versus  nu¬ 
merical  predictions. 

ascribed  to  the  finite  dimensions  of  pressure  taps  and 
of  the  thermocouples.  In  addition  to  the  discrepancy 
in  the  peak  values,  a  second  disagreement  between  ex¬ 
periments  and  numerical  results  is  present  in  the  range 
between  —30°  and  —50°.  However,  numerical  simula¬ 
tions  performed  by  other  authors  for  the  same  test  case 
are  in  very  good  agreement  with  the  present  numeri¬ 
cal  data.  One  possible  explanation  for  the  difference 
with  the  experiment  could  be  attributed  to  some  un¬ 
steadiness  in  the  incoming  shock  position  during  the 
experimental  tests. 

The  numerical  simulation  of  the  experimental  con¬ 
ditions  of  RUN  43  were  performed  using  three  different 
grids.  The  coarser  grid  is  the  one  used  for  investigating 
the  effect  of  the  incoming  shock  position  and  is  com¬ 
posed  of  150  points  in  the  normal  direction,  with  grid 
stretching  close  to  the  wall.  In  the  tangential  direction, 
the  grid  spacing  is  1  degree  in  the  range  [—40°  :  +40°]. 
Then,  the  tangential  grid  spacing  increases  linearly 
up  to  a  maximum  value  of  1.5  degrees  at  the  up¬ 
per  and  lower  boundaries,  so  that  the  grid  dimension 
is  150x143  points.  A  second  mesh  still  contains  150 
points  in  the  normal  direction,  but  the  tangential  spac- 


a)  Mach  number  contours. 


b)  Surface  heat  transfer  distribution. 

Fig.  13  Mach  number  and  logio  of  pressure  for 
RUN  43. 

ing  in  the  range  [—40°  :  +40°]  is  now  0.5  degrees,  with 
maximum  amplitude  of  1  degree  at  the  ends  of  the 
computational  domain  (150x263  points).  For  the  finer 
grid,  300  points  (with  stretching)  where  used  in  the 
normal  direction,  while  the  tangential  spacing  is  0.25 
degrees  between  [—40°  :  +40°]  and  the  maximum  tan¬ 
gential  spacing  is  0.5  degrees  (300x528  points). 

RUN  43  belongs  to  the  family  of  the  Edney-type  IV 
configurations  also.  The  Mach  number  and  logio  of 
pressure  contours  shown  in  figures  13a  and  13b  clearly 
reveal  the  presence  of  a  supersonic  jet  directed  towards 
the  cylinder  wall  and  of  a  normal  shock  that  deceler¬ 
ates  the  flow  before  it  reaches  the  body. 

A  comparison  with  the  pressure  and  the  heat  trans¬ 
fer  at  the  wall  measured  in  the  experiments  is  displayed 
in  figures  14a  and  14b.  Again  the  effect  of  the  grid  re¬ 
finement  is  such  that  a  finer  grid  produces  a  higher 


peak.  The  comparison  with  the  experiments  is  not 
satisfactory  as  far  as  the  peak  values  are  concerned. 
Also  in  this  case,  it  is  possible  to  notice  a  second  dis¬ 
crepancy  in  the  range  [—20°  :  —50°].  The  finer  grid 
shows  the  trend  to  get  closer  to  the  experimental  data 
in  that  area,  arising  the  question  whether  a  further 
grid  refinement  would  be  necessary  to  capture  a  local 
physical  phenomenon.  However,  the  question  about 
the  steadiness  of  the  experiment  remains  open. 

As  mentioned  in  a  previous  section,  in  the  exper¬ 
iments  conducted  at  CALSPAN  it  was  preferred  to 
use  models  with  non-conducting  surfaces.  Then,  in 
those  cases  where  the  resulting  wall  temperature  rises 
were  important,  the  surface  temperature  distribution 
was  tabulated  together  with  the  heat  transfer  rate, 


a)  Surface  pressure  distribution. 


b)  Surface  heat  transfer  distribution. 

Fig.  14  Pressure  and  heat  flux  distribution  at  the 
wall  for  RUN  43.  Experimental  values  versus  nu¬ 
merical  predictions. 
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a)  Surface  pressure  distribution. 


b)  Surface  heat  transfer  distribution. 

Fig.  15  Pressure  and  heat  flux  distribution  at  the 
wall  for  RUN  43  using  a  constant  wall  temperature 
of  300  K  or  the  wall  temperature  measure  in  the 
experiments. 

so  that  it  could  be  used  in  future  numerical  predic¬ 
tions  as  a  boundary  condition.  In  order  to  repro¬ 
ducing  at  the  best  the  experimental  conditions,  the 
numerical  simulations  were  conducted  using  the  wall 
temperature  measured  in  the  experiments  as  a  bound¬ 
ary  condition.  However,  a  comparison  between  two 
computations  performed  using  the  coarser  grid  with 
and  without  the  variable  wall  temperature  condition 
showed  that  the  effect  is  negligible,  as  shown  in  figure 
15. 


for  different  impinging  shock  positions.  Nevertheless, 
computations  are  very  sensitive  to  the  grid  refinement 
and  it  is  very  difficult  to  obtain  grid  convergence  in  the 
area  where  the  wall  pressure  and  the  heat  flux  reach 
their  maximum  values.  In  most  coditions,  the  peak 
levels  measured  in  the  experiments  are  smaller  with 
respect  to  the  numerical  predictions.  This  is  a  very 
critical  point,  since  the  estimation  of  the  highest  heat 
flux  is  fundamental  for  the  design  of  hypersonic  ve¬ 
hicles.  Further  accurate  experimental  and  numerical 
investigations  should  be  necessary  in  order  to  clarify 
such  an  issue. 
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CONCLUSIONS 

Computational  Fluid  Dynamics  is  a  very  useful  tool 
for  simulating  shock-shock  interactions,  since  it  con¬ 
sents  to  predict  the  heat  loads  on  the  body  surface 
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